Quantum coherence and entanglement in the avian compass 
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Tremendous efforts are underway to build technologies that harness the deep quantum phenomena 
of superposition and entanglement. These properties have proven fragile, often decaying rapidly 
unless cryogenic temperatures are used. Could life have evolved to exploit such phenomena [l]? 
Certain migratory birds have the ability to sense very subtle variations in the Earth's magnetic 
field Here we use recent experimental observations ^ together with the well developed 'radical 
pair' model of the avian compass [4], and employ a master equation with various decoherence 
operators in order to examine the system's vulnerability to environmental noise. Remarkably, the 
room temperature noise tolerance in this natural system appears greater than that of the best man- 
made molecular radical ^ or solid state singlet /triplet devices [6j. We find that entanglement, 
though probably not an essential feature of this process, appears to persist to tens of microseconds, 
or more. 



Recently several authors have raised the intriguing pos- 
sibility that living systems may use non-trivial quantum 
effects to optimise some tasks. Studies range from the 
role of quantum physics in photosynthesis [3 [3 |9l fTO] 
and in natural selection itself [11 , through to the obser- 
vation that 'warm and wet' living systems can embody 
entanglement given a suitable cyclic driving [12 . In this 
letter, we examine a particularly important form of nat- 
ural information processing known as magnetoreception 
- the ability to sense characteristics of the surrounding 
magnetic field. There are a several mechanisms by which 
this sense may operate [2]. In certain species (including 
some birds [4 , fruit flies [131 [28 and even plants [14 ), the 
evidence supports a Radical Pair (RP) mechanism, which 
relies on the quantum evolution of a spatially-separated 
interacting pair of electron spins. This is supported by 
results from the field of spin chemistry [T^ [161 (13 UHl [H] 
and recent experiments which were able to demonstrate 
a chemical compass |20j. 

By manipulating a captive bird's magnetic envi- 
ronment and recording its response, one can make 
inferences about the mechanism of the magnetic sen- 
sor [31 [2TJ [22l [23] . Specifically, European Robins are only 
sensitive to the inclination and not the polarization of the 
magnetic field [21 , and this sensor is evidently activated 
by photons entering the bird's eye j22J. Importantly for 
the present analysis, a very small oscillating magnetic 
field can disrupt the bird's ability to orientate ^ [23] . 
It is also significant that birds are able to 'train' to 
different field strengths, suggesting that the navigation 
sense is robust, and unlikely to depend on very special 
values for the parameter in the model 
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All these experiments can be explained with the 
commonly-accepted Radical Pair (RP) model, see Fig.l. 
Here we employ the simplest RP model, considering the 
spins of two electrons [4, 24 and one nucleus of the 
molecule. Absorption of a photon and subsequent trans- 
fer of one electron to an acceptor part of the molecule, 
gives rise to a radical pair initially in the singlet state. 
Due to the spatial separation it now becomes meaning- 
ful to talk about electron spin entanglement. Because of 
the nuclear interaction with one of the electron spins, the 
singlet state is no longer an eigenstate of this Hamilto- 
nian leading to an angle dependent singlet-triplet oscil- 
lation. Recombination occurs either from the singlet or 
triplet state, leading to different chemical endproducts. 
The concentration of those products constitutes a macro- 
scopic chemical signal, correlated to the orientation of the 




bird's eye retina 



FIG. 1: Schematic of the bird's eye. The back of the eye con- 
tains numerous molecules 30 - , fixed with specific orientations 
but not necessarily ordered. In the simplest RP model, each 
such molecule involves three crucial components (see inset): 
there are two electrons, initially photo-excited to a singlet 
state, and a nuclear spin that couples to one of the electrons. 
This coupling is anisotropic, so that the molecule has a direc- 
tionality to it. 
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molecule with the magnetic field. 

We employ the Hamiltonian corresponding to the sys- 
tem once the two electrons have become separated. The 
anisotropic hyperfine tensor coupling the nucleus and 
electron 1, is conveniently written in its diagonal basis 
A = didig{Ax^ Ay^ Az)^ and we assume an axially sym- 
metric (or cigar-shaped) molecule with A^ = 10~^ meV 
and Ax = Ay = Az/2. This is the simplest assumption 
that can provide us with directionality, and we have cho- 
sen the general shape and magnitude of the tensor to be 
consistent with [W. The Hamiltonian is 

i7 = /- A.^i+7B-(^i+^2), 

where / is the nuclear spin operator, Si = {ctx^ Cy^ az)i 
are the electron spin operators {i = 1,2), B is the mag- 
netic field vector and 7 = ^/j^o9 the gyromagnetic ratio 
with /io being Bohr's magneton and g = 2 the g-factor. 
The factor 1/2 in the gyromagnetic ratio accounts for the 
fact that we have a spin one-half system, but we will use 
Pauli matrices such as = diag{l, —1} etc. Here only 
one electron is coupled to one nucleus, whereas the re- 
mote electron is so weakly interacting that we describe it 
as free. Previous authors have considered the case where 
more than one nucleus couples to the system [H [191 E] , 
In the Supplementary Information we show our basic con- 
clusions are not affected by varying the hyperfine tensor, 
adding a second nuclear spin, or even by replacing the 
term completely with a anisotropic electron g-factor. 
Generally, the magnetic field we employ is 

B = 5o(cos(/:?sini^, sin(/?sin'i^, cosi?) 

+ coscc;t(cos0sin6>, sin^sin^, cos^), (1) 

where Bq = 47 jiT is the Earth's magnetic field in Frank- 
furt [3 , and the angles describe the orientation of mag- 
netic field to the basis of the HF tensor. B^i = 150 nT 
is an additional oscillatory field only applied in our sim- 
ulations where explicitly mentioned. For resonant exci- 
tation with the uncoupled electron spin, hco = 2jBo^ so 
that u = cc;/(27r) = 1.316MHz. 

Because of the axial symmetry of the HF tensor we can 
set (f = and focus on the i9 in the range [0, 7r/2] without 
loss of generality. Furthermore, for the oscillatory field, 
we set = 0. 

We model the dynamics of the system with a quantum 
master equation (ME) approach. We employ operators 
representing the relaxation processes; specifically, we in- 
clude two 'shelving states' which represent the system 
having decayed either from an electron singlet state, or 
from one of the triplet states. Ultimately one of these 
two forms of relaxation will occur. The three spins span 
an 8 dimensional Hilbert space to which we therefore add 
two further levels l^) and |T) for the singlet and triplet 
decay outcomes, respectively. The populations of these 
levels will then correspond the singlet and triplet yield. 
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FIG. 2: Angular dependence of the singlet yield in the 
presence of an oscillatory field. The blue reference curve 
shows the singlet yield obtained in the Earth's magnetic field 
with Bo = 47/iT), which is independent of the decay rate k 
for k < 10^ s~^. For better visibility, the blue curve has been 
shifted upwards by 0.001. The red curves show the singlet 
field when a 150 nT magnetic field oscillating with a frequency 
resonant with the Zeeman splitting of the uncoupled electron 
(1.316 MHz) magnetic field is superimposed perpendicular to 
the direction of the static field. We see that this only has 
an appreciable effect on the singlet yield once k is of order 
10^ s~^. Inset: a European Robin (© David Jordan) 



With the usual definiton of singlet \s) and triplet states 
\ti) in the electronic subspace, while | t) and | |) describ- 
ing the states of the nuclear spin, we define the following 
decay operators: 

Ps,^ =\S){s,U 
Pt„^ =\T){to,U 
Pt^,^ =\T){U,U 
Pt_,^ =\T){t.,U 

and similarly for the 'down' nuclear states. This gives us 
a total of two singlet projectors and six triplet projec- 
tors. For simplicity and because this choice corresponds 
exactly to the expression for singlet yield used in previ- 
ous literature, all eight projectors have the same decay 
rate Tp = k. 

For our model we start from a initial density matrix 
p(0) corresponding to the electrons in a pure singlet state, 
and a completely mixed nuclear state, i.e., 

P(0) = idn \S){S\ = \S, i){s, i I + \S, T)(5, T |. 

The decay to the two shelving levels is then described 
using a standard quantum ME with above decay opera- 
tors which effectively discrimate singlet and triplet decay 
events 

p=-^[H,p].+k(Yl P^pPi - \ {Pi^^P + pPiP) ) (2) 
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FIG. 3: Angular dependence of the singlet yield in 
the presence of noise. A) All curves were obtained using 
A; = 10^. The blue curve provides a reference in the absence 
of noise and the red curves show the singlet yield for different 
noise rates given. It is apparent from the plot that a noise 
rate T > 0.1k has a dramatic effect on the magnitude and 
contrast of the singlet yield. 



Note that Eqn.([2| does not contain any environmental 
noise, though this does not alter the estimates of the 
decay rate k (see Supp. Information). In the previ- 
ous literature on the RP model of the avian compass, 
it has been common to employ a Liouville equation to 
model the dynamics. In fact, a term- by-term compar- 
ison of the evolution of the terms of the density ma- 
trix readily confirms that this former approach and our 
ME are exactly equivalent in the absence of environmen- 
tal noise. Both agree with the singlet yield integral, 
^ = /q^('0~ |Trn(p(t))|?/^~)/ce~^*(it, another commonly 
used quantity from the prior literature, when singlet and 
triplet reaction rates are equal. Specifically, the ultimate 
population of our singlet 'shelf l^*) then corresponds to 
However, when we presently wish to introduce vari- 
ous kinds of noise operator, in particular the pure phase 
variant, then the ME approach provides a more intuitive 
framework. 

We now wish to determine an appropriate choice for 
our parameter k in Eqn. |2] In Ref. 3, the authors re- 
port that a perturbing magnetic field of frequency of 
1.316 MHz (i.e. the resonance frequency of the 'remote' 
electron) can disrupt the avian compass. They note that 
this immediately implies a bound on the decay rate (since 
the field would appear static for sufficiently rapid decay) . 
Here we aim to refine this bound on k by considering the 
oscillating magnetic field strength which suffices to com- 
pletely disorient the bird's compass, i.e. 150 nT. (Indeed, 
even a 15 nT field was reported as being disruptive, but 
to be conservative in our conclusions we take the larger 
value here.) To model this effect, we activate the oscilla- 
tory field component defined in Eqn. [l] and examine the 



singlet yield (i.e. the eventual population of shelving state 
\S)) as a function of the angle between the Earth's field 
and the molecular axis. Consistent with the experimental 
work, we find that there is no effect at such weak fields 
when the oscillatory field is parallel to the Earth's field. 
Therefore for our analysis we set the oscillatory field to 
be perpendicular. The results are shown in Figure [2j 
We conclude that if the oscillating field is to disorient 
the bird, as experiments showed, then the decay rate k 
should be approximately 10^ s~^ or less. For higher val- 
ues of k (shorter timescales for the overall process) there 
is no time for the weak oscillatory field to significantly 
perturb the system; it relaxes before it has suffered any 
effect. Such a value for the decay rate is consistent with 
the long RP lifetimes in certain candidate cryptochrome 
molecules found in migratory birds p5] . 

Taking the value k = 10^ s~^, we are able to move to 
the primary question of interest: how robust this mech- 
anism is against environmental noise. There are several 
reasons for dehorence. For example dipole interaction, 
electron-electron distance fluctuations and other parti- 
cles' spin interaction with the electrons will cause deco- 
herence. We describe generic environmental noise with 
a standard Lindblad ME technique [26 , where Eqn. [2] 
above is extended with a dissipator as follows: 

P= Eqn.? +Er^ (l.pLI-1(lIl,p + pLIl;)\3) 

i ^ ^ 

Noise operators Li are (Jy^ (Jz for each electron 
spin individually (i.e. tensored with identity matrices 
for the nuclear spin and the other electron spin). This 
gives a total of six different noise operators Li and we 
use the same decoherence rate F for all of them. We are 
now in a position to determine the approximate level of 
noise which the compass may suffer, by finding the mag- 
nitude of F for which the angular sensitivity fails. This 
is shown in Fig. [3] Conservatively, we can say that when 
F > /c, the angular sensitivity is highly degraded. This 
is remarkable, since it implies the decoherence time of 
the two-electron compass system is of order 100 fis or 
more [32]! To provide context for this number, we note 
that the best laboratory experiment involving preserva- 
tion of a molecular electron spin state has accomplished 
a decoherence time of 80/i5 [5 . 

It is interesting to ask, what is the significance of entan- 
glement between the spins in the avian compass? Having 
inferred approximate values for the key parameters, we 
can plot an appropriate entanglement measure over the 
course of the process, from the initial singlet generation 
to the eventual decay. The metric we use is negativity: 



N{p) = 




where is the trace norm of the partial transpose 

of the system's density matrix. The transpose is applied 
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FIG. 4: A plot indicating the decline and disappearance of en- 
tanglement in the compass system, given the parameter k and 
the noise severity F defined above. Here the angle between 
the Earth's field and the molecular axis in 7r/4, although the 
behavior at other angles is similar. The entanglement metric 
is negativity as defined in the text. 



to the uncoupled electron, thus performing the natural 
partitioning between the electron, on one side, and the 
coupled electron plus its nucleus, on the other. Fig. |4] 
shows how this negativity evolves under our noise model. 
Clearly, the initial singlet state is maximally entangled. 
Under noise, entanglement falls off at a faster rate than 
the decay of population from the excited state. 

The noise model above is of the most general kind, 
representing a unrestricted degradation of the system's 
state due to environmental interactions. We further have 
found that the compass mechanism is almost immune 
to pure phase noise (typically the most aggressive noise 
in artificial systems). Even starting from a fully de- 
phased state (1 5) (5 1 + |to)(^o|)/2, the compass operates 
well. However, in the Supplementary Information we 
show that if such noise were naturally present at a high 
level in the compass (exceeding the generic noise level 
F by more than an order of magnitude) then it would 
render the bird immune to the weak oscillatory magnetic 
fields of Ref. [3 . Thus the sensitivity to oscillatory fields 
implies that both amplitude and phase, and thus entan- 
glement, are indeed protected within the avian compass 
on timescales exceeding tens of microseconds. It is not 
clear why such remarkable protection occurs, but given 
the widely-accepted RP model together with the recent 
experimental data |3j, this conclusion follows. 
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We note that in the final stages of preparation of this 
manuscript, a related work has appeared which consid- 
ers a chemical magnetometer in the context of quantum 
control [27 . 
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THE TERMINOLOGY AND TECHNIQUES USED 
IN THE MANUSCRIPT 

The terminology used in the main paper arises from 
the foundations of quantum mechanics and has more re- 
cently become important in the field of quantum infor- 
mation (QI). Thus we describe the avian compass and 
its dynamics in terms of quantum coherence, environ- 
mental decoherence, and entanglement; moreover we use 
the quantum master equation technique to model the sys- 
tem's dynamics. These terms and techniques are rather 
different to those of the established radical pair (RP) and 
avian compass communities, therefore it is appropriate to 
justify their use and relate them to the established ter- 
minology. 

The primary conclusion of the paper concerns the ex- 
tent to which quantum coherence is protected in the com- 
pass system. By this term, we refer to the purity of the 
system, which is degraded by interactions with the en- 
vironment because these will generally act as a kind of 
measurement of the state (for example, by taking a quan- 
tum of energy from the state and thus leaving it in a spe- 
cific low-energy state). Since the 'outcome' of this mea- 
surement is lost to the environment, the system's state 
becomes a mixture of the multiple possible states - gener- 
ically this loss of purity is called decoherence. 

In the applied QI field, much effort is dedicated to the 
preservation of quantum coherence. In order for quan- 
tum coherence to be maintained, it is essential to prevent 
all kinds of interaction, including those that would flip 
spins, allowing them to relax and lose energy, and also 
those interactions which would merely alter the phase re- 
lationships in the state (dephasing noise). Thus while the 
concept of spin relaxation is a typical signature of general 
coupling to the environment, we must also examine the 
effect of pure phase noise; the latter alone can suffice to 
completely degrade a pure state such as the singlet \s) 
into a completely incoherent mixture of \s) and \to) and 
would therefore undermine our conclusion. In order to 
examine this point, and establish whether full coherence 
is indeed maintained, it is necessary to construct a equa- 
tion for the dynamics of the system where pure phase 
noise alone occurs. This is most naturally done with 
the quantum master equation (QME) formalism and the 
corresponding Lindblad noise operators - however, it is 
important to note (as mentioned in the main paper) that 
the QME is actually identical to the more conventional 
Liouville equation in the case of zero environmental noise. 
We select the QME as the most suitable tool for examin- 
ing specific decoherence models, including both general 
spin relaxation and those that specifically cause only de- 
phasing. 
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Finally we note that, because the RP involves two elec- 
trons, and these are spatially separated, it is meaningful 
to use the quantity 'entanglement' as one specific char- 
acterization of their mutual state. The established avian 
compass literature would tend to refer to the related con- 
cept of spin correlation^ however, entanglement is subtly 
and importantly different in that it captures precisely the 
non-classical correlations between spins - it is impossible, 
for example, to create entanglement between two sites 
simply by exchanging classical information as to how the 
spins should be prepared. We chose to plot the entan- 
glement in the avian compass because it is a key quan- 
tity (arguably, the essential quantity) in the quantum 
mechanics, and therefore a reader from that community 
will naturally wonder about this aspect. However, we 
find that while entanglement does persist for a remark- 
ably long time, this is likely to be merely an implication 
of the long coherence time rather than a key property of 
the compass mechanism. 
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FIG. 5: Angular dependence of the singlet yield at /c = 
lO'^ s"^ in the presence of the oscillatory field for different 
pure dephasing rates F^. This is to be compared with the 
k — 10^ s~^ line in Fig. 2 of the main paper. See text for an 
explanation. 



SENSITIVITY TO OSCILLATORY FIELD IN THE 
PRESENCE OF PURE DEPHASING 

Interestingly, if we begin the simulation with a com- 
pletely dephased state: + |to)(to|)/2, the classical 
correlations are still sufficient for achieving adequate an- 
gular visibility and neither quantum phase coherence nor 
entanglement seems to be a prerequisite for the efficiency 
of the avian compass. 

To explore this idea further, we would like to study 
'pure dephasing' occurring during the singlet-triplet in- 
terconversion. In essence we use energy conserving noise 
operators, Eqn. Q, which are known to be the domi- 
nant source of decoherence in so many other artificially 
made quantum systems. By applying this specific noise, 
we confirm that the compass mechanism's performance 
is essentially immune, while of course the coherence of 
the quantum state of the electrons would be degraded. 

One might be inclined to conclude that, if pure de- 
phasing noise is indeed dominant, then the avian com- 
pass need not protect quantum coherence for the long 
time scales suggested in the main paper. But crucially, 
we also show that if such noise were naturally present at 
a high level in the compass (exceeding the generic noise 
level F by more than an order of magnitude) then it would 
render the bird immune to the weak oscillatory magnetic 
fields studied by Ritz et al. [1 . Thus the sensitivity to 
oscillatory fields implies that both amplitude and phase, 
and thus entanglement, are indeed protected within the 
avian compass on timescales exceeding tens of microsec- 
onds. 

Since the electron spin singlet state is not an eigen- 
state of the Hamiltonian, the dephasing operators will 
be different from the ones mixing the phase of the singlet 
and triplet state within the electronic subspace. Instead, 



we replace the previously defined noise operators of Li of 
Eq. (3) by appropriate dephasing operators as follows: 
we treat the remote electron and the eletron nuclear spin 
subsystem separately. Within both subsystems, we de- 
fine dephasing operators 

= -l(/4-2|A,)(A,|), (4) 

where {|Ai)} are the set of normalised eigenvectors of this 
subsystem. This results in two dephasing operators for 
the remote electron (these can be combined to a single 
(Jz operator rotated with the field) and four operators 
for the eletron nuclear spin subsystem. Each of these 
dephasing operators corresponds to fluctuations of one 
of the (subsystem's) energy levels. 

Strikingly, the singlet yield is entirely unaffected by 
this particular kind of noise, i.e. it is entirely independent 
of the dephasing rate F^. Thus, a curve obtained with 
this model coincides perfectly with the reference curve 
of Fig. 3 of the main paper. However, we show in the 
following that the dephasing rate of this model can be 
at most ten times faster than the generic noise rate to 
retain sensitivity to the oscillatory field. 

Fig. [5] shows the singlet yield as a function of 9 for dif- 
ferent pure dephasing rates F^. Pure phase noise would 
actually protect the compass from the harmful effect of an 
applied oscillatory field (by suppressing the Rabi oscilla- 
tions caused by such a field). We see that an aggressive 
pure dephasing rate of 1/Tz = 10 /is almost completely 
recovers the reference curve (corresponding to a noise- 
free system without oscillatory field). 
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GENERALITY OF THE MODEL 

European robins can adjust to different strengths of the 
magnetic field. This impUes, that their chemical compass 
cannot depend on the absolute value of the the electron 
spins hyperfine (HF) coupling. We here argue that the 
conclusion presented in the main paper is not only robust 
against variations in the field strength, but also does not 
depend on the specifics of the anisotropic interaction be- 
tween electron 1 and its local environment. (Of the two 
electrons constituing the RP, only one of them can be 
significantly coupled to nuclear spin(s) because the res- 
onance effect [l] occurs at exactly the frequency of the 
Zeeman splitting of an unperturbed electron.) 



'Disc-shaped hyperfine tensor' 

We have explicitly checked that our conclusions do not 
rely on the specifics of the HF tensor given in the main 
text. Specifically neither a weaker or stronger coupling, 
nor a different HF tensor symmetry changes our obser- 
vations about the shortest time scale required for the 
process to be sensitive to the oscillatory field. Neither do 
these factors change the maximal tolerable environmen- 
tal noise rates to maintain a pronounced signal contrast. 

Here we describe one of the variants we have studied: a 
'disc-shaped' HF-coupling tensor with coupling strengths 
that are different from the ones presented in the main 
paper. The largest coupling here is half as strong as the 
one presented in the text: = 0.5 x 10~^ meV, Ay = 

Fig. |6] shows that we obtain exactly the same quali- 
tative behaviour as for the parameters used in the main 
text, the only difference being in the different shape of 
the singlet yield curve (as expected for a different geom- 
etry of HF tensor). The oscillatory field does not have 
an infiuence on the singlet yield for k > 10^ s~^, while a 
noise rate of F > 0.1/c leads to a dramatic reduction of 
contrast. 



Anisotropic g-factor 

Instead of coupling to nuclear spin(s), an anisotropic g- 
factor for one of the electrons could also lead to a singlet 
yield with an angular dependence [2 . We here show that 
our main conclusion remains valid under this different 
Hamiltonian of the RP. 

For illustration purposes, we consider a prounounced 
anisotropy: the g-factor in z-direction is 0.8 x 2, i.e. 80% 
of the free electronic g-factor of ^ = 2, while in x and 
^/-direction we have 0.3 x 2. A smaller anisotropy works 
as well, but gives less overall signal contrast (even in 
the unperturbed) scenario. We present results for these 



parameters in Fig. [7| confirming that the qualitative 
behaviour is still the same under these assumptions. 



Markovian noise model 

The noise model we employ is Markovian, implying 
that the information lost to the environment will not re- 
turn to the system (as it would, wholly or partially, if 
the environment was entirely constituted by other iso- 
lated qubits that could effectively act as a memory). This 
Markovian assumption is the natural choice to make for a 
room temperature system whose detailed structure and 
environmental interactions are unknown. One can im- 
mediately think of 'likely suspects' for decoherence in 
the present system that would be Markovian: for exam- 
ple, modulation in the electron-electron separation due 
to vibrations in the supporting matrix. But more impor- 
tantly, while one can reasonably speculate that there may 
also be non-Markovian processes occuring, between the 
core system and other degrees of freedom that are well 
isolated from the true 'bath' degrees of freedom, this does 
not undermine our result. In essence, one would 'draw a 
dotted line' around the larger system, comprised of the 
core electron pair together with the non-Markovian envi- 
ronment, calling that the avian compass (although noting 
that only the core part plays an active role). One would 
conclude that this larger system is, again, remarkably 
well insulated from the environment. 



RP PAIR MODEL WITH 2 NUCLEAR SPINS 

In the model described in the main paper, there is a 
single nuclear spin coupled to one of the electrons. How- 
ever, previous publications have studied the case where 
more than one nuclear spin is present [U El [4]. There- 
fore, it is interesting to check whether the addition of a 
nuclear spin will alter our conclusions. The Hamiltonian 
now gains an additional coupling term, 

i7 = A • Ai • ^1 + /2 • A2 • ^1 + 7B • (^1 + ^2), 

where Ai is the HF tensor of the main text and A2 = 
2/3 Ai. Here we choose a second rugby shaped HF ten- 
sor oriented parallel to the first. Note that we have also 
considered different relative coupling strengths and ge- 
ometries (such as a pancake shaped and rugby shaped 
tensor), but again, these choices do not infiuence our core 
conclusions. Results for the particular choice of parame- 
ters described above are shown in Fig. [S] 
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FIG. 6: Results for a disc-shaped HF tensor. The graph on the left is precisely analogous to Fig. 2 of the main text, showing 
angular dependence of the singlet yield in the presence of the oscillatory field for different decay rates k. The graph on the 
right corresponds to Fig. 3 of the main paper, showing angular dependence of the singlet yield in the presence of environmental 
noise [obtained with Eq. ^] for different noise rates F. Although the curve shapes differ, the levels of contrast are the almost 
the same as the 'cigar' shaped HF model in the main paper. 
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FIG. 7: Results for a model in which an anisotropic g- factor replaces the role of the nucleus to break the RP symmetry. The 
graph on the left is precisely analogous to Fig. 2 of the main text, showing angular dependence of the singlet yield in the 
presence of the oscillatory field for different decay rates k. The graph on the right corresponds to Fig. 3 of the main paper, 
showing angular dependence of the singlet yield in the presence of environmental noise for different noise rates F. Although 
the curve shapes differ, remarkably the levels of contrast remain similar to those of the conventional model in the main paper. 



OSCILLATORY FIELD SENSITIVITY IN THE 
PRESENCE OF NOISE 



In the first stage of the argument in the main paper, 
we show how the bird's sensitivity to field angle is de- 
graded by an applied oscillatory field (Fig. 2, main pa- 
per). There we set the environmental noise operators to 
zero - this is of course unrealistic, and one might worry 
that finite noise could ultimately lead to less dramatic 
conclusions regarding protection of quantum coherence. 
However, here we show that the reverse is true: the addi- 



tion of finite noise leads to a suppression of the harmful 
effect of the oscillatory field, and thus would lead one to 
infer a still more dramatic coherence time. Therefore, by 
setting the noise to zero we make the more conservative 
assumption. 



We confirm this by employing the following master 
equation (which is exactly Eq. (3) of the main paper). 
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FIG. 8: Results for two nuclear spins coupled to electron 1. The graph on the left is precisely analogous to Fig. 2 of the main 
text, showing angular dependence of the singlet yield in the presence of the oscillatory field for different decay rates k. The 
graph on the right corresponds to Fig. 3 of the main paper, showing angular dependence of the singlet yield in the presence 
of environmental noise for different noise rates F. Although the curve shapes differ somewhat, nevertheless the parameters 
corresponding to serve degradation remain the same as those in the main paper (k = 10^ s-\ F = 100 /is). 



while at the same time applying the oscillatory field: 

p = -'-[H,p] + k(j2 PipPi - \ (p\Pip + pplp^ 

6 , ^ 

+ {upL\ - - (L\Lip + pL\L,) 



(5) 



The left panel of Fig. |9] corresponds to Fig. 2 of the 
main paper with such additionally applied noise at a rate 
of F = 0.1/c. Note that the noise rate is adjusted with 
/c, as we wish to concentrate on the combined effect of 
the oscillatory field and 'mild noise' over the timescale 
of the process. We infer that the singlet yield contrast 
is only significantly reduced when k < 10^ s~^, meaning 
the addition of noise to Fig. 2 in the main paper does 
not influence our conclusion. 

To establish that this conclusion is robust with respect 
to the particular choice of F = 0.1/c, we now investigate 
different noise rates for a fixed value of k = 10^ s~^. 
In the right panel of Fig. |9j we show a pair of curves 
for each noise rate: one obtained with and one without 
the oscillatory field. We make two observations: firstly, 
the two curves of each pair nearly coincide, indicating 
that there is only a small sensitivity to the oscillatory 
field indepedent of the noise rate. Secondly, the sensitiv- 
ity decreases even more when the system is subjected to 
more aggressive noise. 

In fact, this second observation is not surprising: we 
would not expect the additional noise terms to lead to 
a greater sensitivity to the oscillatory field. The singlet 
yield signal is corrupted by spin flips of the remote elec- 
tron in the presence of the oscillatory field, which essen- 
tially causes the spin to Rabi flop, thus requiring quan- 



tum coherence. Therefore, any noise acting on the remote 
electron will decohere the spin and thus suppress rather 
than enhance the effect of the oscillatory field. This can 
also be seen in the right panel of Fig. [9j the red curves 
are very similar but not quite identical, the agreement 
is better for the green curves and almost perfect for the 
blue one. 
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FIG. 9: Left: singlet yield as a function of in the presence of the resonant oscillatory field with additionally applied noise at 
a rate F = 0.1/c (see text). Right: The singlet yield as a function of 0^ this time with k fixed at A; = 10^ s~^ and for different 
noise rates. For each noise rate, there are two curves, one with and one without the oscillatory field. However, as explained in 
the text, these can only be properly distinguished for the case of weak noise (red curves). 



